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Measurements  of  the  sound  radiation  patterns  in  water  from 
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tions made  in  the  theory,  it  was  found  that 'the  major  features 
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I.   INTRODUCTION 

A  high  frequency  broad-band  transducer  which  has  a  broad 
radiation  pattern  in  a  half-space  is  required  for  acoustic 
tracking  of  underwater  vehicles  in  situations  where  the  vehicle 
carries  the  sound  source  and  the  receivers  are  fixed.  [1]  The 
underwater  ranges  at  the  Naval  Undersea  Warfare  Engineering 
Station  (NUWES) ,  Keyport,  Washington,  are  of  this  type.   During 
tests  of  production  torpedoes  and  new  torpedo  designs,  the  tor- 
pedo is  tracked  in  three  dimensions  as  it  moves  through  the 
water . 

The  tracking  is  accomplished  using  an  array  of  hydrophones 
mounted  on  the  bottom  of  the  acoustic  range  which  receives 
acoustic  pulses  radiated  periodically  from  a  projector  mounted 
on  each  vehicle  on  the  range.   Currently,  an  acoustic  projector 
mounted  flush  with  the  bottom  of  the  torpedo  sends  out  a  7  5kHz 
coded  pulse  that  is  picked  up  by  the  hydrophones.   A  shore- 
based  computer  receives  the  information  from  the  range  hydro- 
phones and  computes  the  location  and  speed  of  the  torpedo 
during  the  run. 

For  the  NUWES  Keyport  ranges,  the  signal  to  noise  ratio  re- 
quirements and  the  expected  maximum  transmission  loss  between 
source  and  the  receivers  places  a  minimum  requirement  on 
source  level  which  is  quite  high.   In  addition,  the  pitch  and 
roll  motions  of  the  torpedo  and  the  effect  of  surface  reflec- 
tion of  the  sound  has  led  to  a  specification  of  a  broad  beam 
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width.   Specifically,  current  requirements  are  for  a  uniform 
source  level  within  a  conical  volume  about  150°  wide  centered 
on  the  vertical  axis  of  the  transducer. 

The  projectors  are  required  to  operate  acoustically  un- 
loaded (in  air)  95  percent  of  the  time,  since  the  unit  is 
turned  on  before  launch  and  continues  running  after  recovery. 
Another  requirement  placed  on  the  projector  to  be  used  in 
higher  speed  vehicles  is  that  the  radiating  surface  be  flush 
with  the  body  of  the  vehicle. 

These  various  requirements  place  extreme  demands  on  the 
physical  properties  of  materials  used  in  the  transducer  and 
create  difficult  problems  of  building  a  transducer  which  also 
has  reasonable  life  expectancy  and  reasonable  cost. 

Another  approach  to  this  problem  is  the  use  of  a  flexural 
wave  transducer  as  a  broad-beam,  broad-band-width  source  which 
would  be  readily  adaptable  to  flush-mounting,  a  necessary  re- 
quirement for  high  speed  torpedoes.   The  possibility  of  using 
the  torpedo  hull  itself  as  the  radiating  surface  is  especially 
attractive  for  the  low  drag,  laminar  flow  shapes  expected  to 
become  used  in  the  future.   Present  designs  use  a  transducer 
which  is  inserted  into  a  recess  in  the  hull.   Although  this  is 
relatively  smooth,  it  is  not  adequately  smooth  for  ultra-low 
drag  laminar  flow. 

Some  early  work  in  this  area  was  done  by  Barlow  at  NUWES 
(private  communication) .   The  first  previous  work  at  the  Naval 
Postgraduate  School  on  this  subject  was  done  by  LT  Sevdik  [1] . 
His  calculations  were  based  on  classical  plate  theory  as  de- 
scribed in  Morse  and  Ingard  [2]  and  Malechi  [3] ,  which  does 
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not  take  into  account  of  effects  of  transverse  shear  and 
rotatory  inertia.   Among  the  results  of  Sevdik's  work  was  the 
noting  of  damping  of  the  plate  natural  frequencies  by  fluid 
loading  and  of  inconsistencies  in  the  change  of  flexural  wave 
phase  speed  with  frequency  for  transducers  of  differing 
thickness. 

Mindlin  [4]  includes  these  effects  in  a  three-dimensional 
plate  theory  and  shows  that  classical  plate  theory  begins  to 
differ  markedly  from  three-dimensional  theory  at  frequencies 
above  that  at  which  the  flexural  wavelength  is  about  ten  times 
the  plate  thickness.   His  theory  gives  f lexural-wave  speed 
values  for  a  given  ratio  of  plate  thickness  to  flexural  wave 
length  which  are  lower  than  those  from  classical  plate  theory. 

Feit  [5]  considers  the  case  of  an  infinite  plate  excited 
by  a  simple  harmonic  point  source  and  which  is  fluid  loaded  on 
one  side.   He  discusses  the  significant  changes  in  the  theoreti- 
cal radiation  patterns  that  result  at  frequencies  above  and 
below  the  critical  frequency,  that  frequency  at  which  the 
flexural  wave  speed  in  the  plate  equals  the  speed  of  sound 
in  the  fluid.   His  work  is  based  on  the  Mindlin-Timoshenko 
plate  theory. 

An  experimental  investigation  was  conducted  of  the  effects 
of  fluid  loading  on  axially  symmetric  flexural  waves  in  a 
circular  aluminum  plate  25  in.  in  diameter  and  5/16  in.  thick 
by  LT  Jarvis  [7]  in  1977.   Jarvis '  results  are  in  good  agree- 
ment with  the  Mindlin-Timoshenko  plate  theory  for  the  speed 
of  the  flexural  wave  as  a  function  of  frequency.   Jarvis'  re- 
sults clearly  show  a  change  in  the  radiation  resistance  as 
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frequency  increases,  as  evidenced  by  the  increased  damping  of 
the  plate  natural  frequencies  with  increasing  frequency.   At 
frequencies  below  the  critical  frequency,  the  fluid  adds  iner- 
tia which  causes  a  lower  flexural  wave  speed.   Above  the  critical 
frequency,  the  wave  speed  increases  abruptly  to  about  the  value 
for  the  unloaded  plate,  indicating  that  the  radiation  reactance 
(inertial)  becomes  very  much  smaller  above  the  critical  fre- 
quency . 

Until  the  work  of  Alper  and  Magrab  [6] ,  a  "closed  form" 
solution  of  the  coupled  vibration  of  a  finite  circular  plate 
clamped  in  an  infinite,  rigid  baffle  had  not  been  obtained  for 
either  classical  plate  theory  or  the  Mindlin-Timeshenko  plate 
theory. 

The  first  treatment  of  this  problem  was  given  by  Lord  Ray- 
leigh  [8]   Mindlin  [10]  andDeresciewicz  [9]  have  shown  that  the 
normal  modes  and  natural  frequencies  of  a  plate  of  finite 
thickness  are  different  and  more  numerous  in  a  given  frequency 
band  that  are  given  by  classical  plate  theory. 

Reisman  and  Greene  [11]  have  studied  the  "in  vacuo"  response 
of  clamped  circular  plates  to  suddenly  applied  loads  using  the 
Mindlin-Timoshenko  plate  theory.   They  found  that  except  for 
the  first  three  or  four  natural  frequencies  of  very  thin 
plates  (diameter  to  thickness  ratio  40  or  more)  there  was  no 
similarity  between  the  natural  frequency  spectra  of  the  Mindlin- 
Timoshenko  theory  and  classical  theory.   The  difference  in 
transient  responses  were  also  significant. 

In  the  Alper-Magrab  paper  the  acoustic  pressure  distribu- 
tion is  obtained  in  a  half  space  bounded  by  an  infinite  rigid 

14 


baffle  containing  a  clamped  circular  plate  which  is  driven 
from  the  "in  vacuo"  side  by  a  simple  harmonic  force  having  an 
arbitrary  spatial  distribution.   A  closed  form  solution  is 
obtained  which  is  valid  for  both  the  near  and  far  fields.   The 
analysis  uses  oblate  spheroidal  coordinates  to  obtain  the  solu- 
tion of  the  wave  equation  in  the  fluid.   These  coordinates  are 
chosen  because  a  bounding  surface  of  one  of  the  coordinates  is 
represented  entirely  by  the  circular  plate,  and  a  bounding  sur- 
face of  another  coordinate  is  represented  entirely  by  the  in- 
finite baffle  with  a  circular  opening  the  same  size  as  the  plate 
diameter.   The  third  coordinate  is  identical  to  one  of  the 
polar  coordinates,  the  coordinates  in  which  the  plate  motion 
is  obtained.   As  a  result,  the  otherwise  fairly  complicated 
coupled  problem  is  reduced  to  two  simpler  boundary  value  pro- 
blems which  can  be  solved  separately  and  the  solutions  coupled 
by  the  continuity  condition  at  their  common  boundary.   The  re- 
sulting solution  for  the  acoustic  field  is  given  by  an  expansion 
in  spheroidal  wave  functions.   Pressure  at  any  point  in  the 
field  can  be  evaluated  by  inserting  the  value  of  the  wave  func- 
tions at  the  point. 

This  is  in  contrast  to  previous  known  attempts  to  solve 
this  problem,  all  of  which  have  resulted  in  the  need  to  perform 
difficult  double  (surface)  integrations  numerically  with  each 
integration  valid  for  only  one  point  in  the  field.   Furthermore, 
all  known  previous  work  has  been  restricted  to  the  classical 
plate  theory. 

In  the  following  pages,  a  brief  review  is  presented  of 
Alper's  theoretical  work  as  reported  in  his  doctoral  thesis  [16] 

15 


and  from  the  point  of  view  of  using  the  results  as  a  design 
tool  for  a  transducer. 

As  part  of  the  present  thesis,  experimental  work  was  conducted 
as  a  test  of  the  Alper-Magrab  theory,  using  measurements  of  the 
radiated  beam  pattern  from  a  transducer.   These  results  are 
presented  and  compared  with  Alper ' s  calculation. 

Finally,  a  modification  of  the  Alper-Magrab  computer  pro- 
gram is  included  as  an  appendix  in  a  form  suitable  for  running 
on  the  Navy  Fleet  Numerical  Weather  Facility  Computer  (CDC- 
6500-"HAL")  . 
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II.  '  THEORETICAL  CONS  IDE-RAT  IONS 

A.   APPROACH 

An  outline  of  the  approach  to  the  problem  [16]  of  calculat- 
ing the  sound  field  from  the  flexual  vibration  in  a  clamped 
edge  plate  is : 

(1)  The  wave  equation  governing  the  motion  in  the  acoustic 
fluid  is  transformed  into  oblate  spheroidal  coordinates. 

(2)  Using  the  Mindlin-Timoshenko  plate  theory,  the  eigen 
frequencies  and  eigen  functions  of  the  plate  are  determined 

for  a  clamped  circular  plate  using  polar  cylindrical  coordinates 

(3)  The  solution  to  the  forced  vibration  problem  is  ex- 
pressed as  the  collection  of  all  possible  combinations  of  the 
plate  eigenfunctions  and  is  introduced  into  the  appropriate 
nonhomogeneous  differential  equation  governing  the  plate  motion 
which  contains  a  term  representing  the  surface  force  acting  on 
the  plate  [the  surface  force  is  the  sum  of  the  applied  force 
and  the  fluid  pressure  (in  terms  of  velocity  potential  expressed 
by  the  means  of  spheroidal  wave  functions) ] .   The  resulting 
equations  are  then  reduced  to  algebraic  equations  in  the  two 
sets  of  arbitrary  coefficients  (plate  motion  and  fluid  motion) 
by  using  the  orthogonality  property  of  the  "in  vacuo"  eigen- 
functions of  the  plate. 

(4)  Because  the  velocities  of  the  plate  and  of  the  fluid 
must  equal  each  other  at  the  interface  of  the  plate  on  fluid, 
the  appropriate  derivatives  of  the  plate  displacement  and  fluid 
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velocity  potential  are  equated.   Orthogonality  of  the  spheroidal 
^ave  functions  then  yields  another  set  of  algebraic  equations. 

(5)   Combining  (3)  and  (4)  and  solving  the  result  gives 
equations  which  express  the  pressure  at  every  point  in  terms  of 
spheroidal  wave  functions,  thus  leading  to  a  knowledge  of  pres- 
sure at  the  point  of  interest  at  a  given  frequency. 

B.   SOLUTION 

The  following  steps  are  from  the  work  of  Alper  [16] .   As 
seen  in  Figure  1  a  circular  plate  is  clamped  in  an  infinite  rigid 
baffle,  represented  by  the  x,  y  plane.   The  upper  part  is  filled 
with  fluid  and  there  is  vacuum  in  the  lower  part.   A  spatially 
arbitrary,  harmonic,  transverse  force  is  applied  to  the  "in- 
vacuo"  side  of  the  plate.   The  objective  is  to  calculate  the  far 
field  pressure  in  the  acoustic  fluid  produced  by  the  forced  motion 
of  the  plate.   The  equation  of  motion  in  the  fluid  is 
2V2    =Q2   3  2  0 

where 
**■  =   ka,  frequency  parameter 
^  =  uj"L,  dimensionless  time. 
0  =  dimensionless  velocity  potential 

p   -  Si         £0        dimensionless  pressure        (2) 

V   =   -  Q.  V0         dimensionless  velocity       (3) 
We  will  look  for  harmonic  solutions  of  the  form 

0  =  fi   e  iZ  (4) 

where 
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0    =  magnitude  of  velocity  potential  and  is  a  function  of 

space  coordinates  only. 

Hence,  putting  Eq.  C4J  into  Eq.  (2)  gives 

p  =  ifl?^eiT  =  peit  (5) 

Again,  p  =  i_ft.£)   is  a  function  of  space  coordinates  only. 
Also,  putting  Eq .  (4)  into  Eq.  (1)  gives 

aV(0eir)  -H2(i2$  eir  )  =  0 

(a2V2  +  -0-2)  0   =   0  (6) 

Eq.  (6)  is  the  dimensionless  Helmholtz  equation.   Its  solutions 
must  satisfy  the  boundary  conditions  in  the  fluid.   Except  for 
the  conditions  at  the  plate,  these  are: 

(a)  The  fluid  velocity  normal  to  the  rigid  baffle  must  be 
zero. 

(b)  Only  outgoing  waves  are  accepted. 

Oblate  spheroidal  coordinates,  which  are  shown  in  Figure  2 
[Van  Buren,  King,  Baier,  and  Hanish,  Ref .  17],  are  the  most 
suitable  for  these  boundary  conditions. 

When  the  confocal  ellipses  and  hyperbolas  having  an  inter- 
focal  distance  of  plate  diameter  2a  are  rotated  about  the  minor 
axis  (z  axis) ,  flattened  ellipsoids  and  hyperboloids  of  one 
sheet  are  obtained.   Under  the  transformation  to  dimensionless 
coordinates  V  and  T. 

x  =  a  [(l-^2)  (\2+±)}1/2  cos9 

y=a  [(1-7J2)  (  |2+1)]1/2     sin9  (7) 

z  =  a  [y\]  . 

The  <!=   constant  surface  is  an  ellipsoid  of  resolution.   The 
cross  section  perpendicular  to  the  x  or  y  axis  is  an  ellipse, 
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while  that  perpendicular  to  the  z  axis  is  a  circle.   In  the  xy 
plane  the  eccentricity  of  the  ellipse  is  -1-  and  the  semimajor 
axis  is   2.  •   Hence  the  bigger  £  is,  the  larger  the  ellipse, 
and  the  more  nearly  circular  it  becomes.   The  surface^  =  0  is  a 
circular  disk  of  radius  a  in  the  xy  plane.   The  surface^=o° 
is  a  sphere  of  radius  °o   .   The  1)  =  constant  surface  is  a  hyper- 
boloid  of  revolution  of  one  sheet  with  anasymptotic  cone  whose 
generating  line  passes  through  the  origin  and  is  inclined  at  the 
angle T)=  cosV  to  the  z  axis.   The  cross  section  perpendicular  to 
the  z  axis  is  an  ellipse,  sections  perpendicular  to  the  x  or  y 
axis  are  hyperbolas.   The  surface V=    0  is  the  xy  plane  except 
for  the  circular  disk  £_=  0.   The  surface  1^1  =  1  is  the  z  axis. 
The  intersection  of  V   and  ^  constant  surfaces  is  a  circle.   The 
angle  9,  measured  from  the  x  axis  towards  the  y  axis  determines 
a  point  on  this  circle. 
Domains  are: 

Eq.  (23)  becomes,  in  oblate  spheroidal  coordinates, 

„i  1 

via 


0  =  0 


(8) 

The  general  solution  of  this  by  separation  of  variables  gives 


3  .s  c-^^)Rnc-^i;©.c^) 


o: 
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where 


and 


S   (-in,  Ti  )  and  R   (-iil,  i?)   represent  angular  and  radial 

functions  respectively.   Putting  the  general  solution  Eq.  (9) 
into  the  Helmholtz  equation  Eq.  (8)  leads  to  two  ordinary- 
differential  equations. 


d_ 


^SUa,i) 


+ 


A  +SLV 

mn    '-      \—yj 


$t*n>° 


(10) 


i_ 


r+. 


RJ-'^'l>°  (id 


X 


where  A   is  the  separation  constant  and  m  =  0,1,2,3 

ran  c 

After  mathematical  calculations 


/N 


0) 


(M) 


^  =  ^n^nC-^^)^o(-^'1)©.C9) 


(12) 


(n-m)  are  even  values 


where 


A    =  constant, 
mn 


oo 


SmnU1  ^  V  "  Z-.     S 


mn 


(-iflO  P™   (7)) 
m- 


1  T 


is  the  oblate  angular  function  of  the  first  kind,  and 

is  the  associated  Legendre  function  of  order  m  and  degree  (m+Q) 
From  [ref.  16,  pp.  17-23]. 

R^}  (-iil,  il  )  =  RCi}  C-iil,  iS  )  -  i  R^}  C-iil,  if  )  • 
mn  I  mn  i  mn  c 
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In  the  development  of  Eq.  (.12)  some  of  the  possible  solutions 

are  rejected.   The  coefficient  of  the  second  independent  solu- 

(2 ) 
tion  of  Eq  (4)  ,  which  is  S    (-ill,  7j  )       ,  is  chosen  to  be 

zero,  because  at  V)  =  1  this  function  goes  infinity.   The  coef- 
ficient of  the  other  independent  solution  combination  of  Eq.  (11), 


which  is  R(3)  (-ill,  i  ?  ) 
mn  I 

only  outgoing  waves  are  considered.   R    &  R    (-i.fi,  i^  ) 


,  is  also  chosen  to  be  zero,  because 

"5  ] 


be  calculated  in  power  expansion  given  by  Flammer  [ref .  18] . 
Also  by  expanding  Eq.  (3)  in  oblate  spheroidal  coordinates  yields 

12  ^ 


v= 


1         u 


+«- 


(13) 


where  ©„,  ,   e    and  e    are  the  unit  vectors  normal  to  the 
**!  ,  J    and  Q  surfaces,  respectively.   The  velocity  component  normal 
to  the Y surface  is 


±2£   ii 

The  baffle  plane  corresponds  to  V 
condition  requires  that 


=  0 


(14) 
The  rigid  baffle 


Vv 


=  0 
7=0 


^=o 


Using  Eq.  (12)  for  Eq .  (15)  yields 


(15) 


d  5 

r 


0) 


7&X 


=  0 

7=o 


(16) 
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After  using  the  power  expansion  for  S     (-iIL,V  ) ,  odd 
values  of  (n-m)  simply  could  not  satisfy  Eq.  (16) .   Therefore 
odd  values  of  (n-m)  are  rejected  and  hence  it  is  found  that 
Eq.  (12)  exists  only  for  even  values  of  (n-m) . 


Summation  of  all  acceptable  solutions  yields 

OO         $0 


A 

t 


A    S(1) 

mn  mn 


(-li^R^t-in,;?)  ®  (G)  .      (17) 
L     mn       c    m 


m=0    n=m,m+2 ,m+4 


In  the  case  of  axisymmetric  motion,  m  =  0  and  then  the 
far  field  solution  becomes 


of 


A1 


CO 


0/2      (1) 
'-1'    AonSon   (-^7' 


n=0,2,4 


(18) 


C.   MINDLIN-TIMOSHENKO  PLATE  ANALYSIS  [REFS  6  and  16] : 

The  dimensionless  equations  of  motion  in  polar  coordinates 
for  simple  harmonic  excitation  are: 


K 


(v1  ^  -Ui.kv^r  ($+-&)¥<Ls?& 


=  0 


(19a) 


"K 


=0  d9b: 
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%r+h 


-,i 


Co   SL 


h 


-sMHpaM) 


^ 


(19c) 


and 


x        3u         u.      W/-^9 


1      ,   Z    T. 
where  _TL  =  n.  u. 


or(u,e)= 


dimensionless  transverse  displacement 


vu,9  'Su'-  u2>u   ui^e! 

Y  =  Vavrrr 

k  =  empirical  she 


ar  coefficient  for  plate. 


A 


/\ 


To-  'vand  T&^u/^/are  the  angles  of  rotation  of  the  nor- 
mal to  the  plate's  middle  surface  a  result  of  bending  in  the  u 
and  "  direction  respectively. 

r 


u  = 


a 


9Cu;^/=   spatially  distributed  harmonic  external  load 

"N  A 

PCLCu,&)=  reaction  of  the  fluid  pressure  to  Q   C^y^V 

A        A 

A 

For   Q  =  0  (free  vibration)  solutions  to  Eq.  (19)  are 


m     m(TA*X+(r1-\)&+±-M 
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A  /\  A» 

NAT  =    ^7  +  W^     , 

A  A  A. 

If    WJ"  ,     \jJ~         and  H      satisfy 

(  VU/e  +  (X,  )  \at  -  o 

<  N7U^  +  <-)  *r  -  o 

(     \7^  +  <    )       H  =0  (21) 

where  ,  ( ,  x        v     r  \ 


(   <T,  <T,  )    =  lfLi_lliJ  (21a) 

Solutions  to  Eq.  (21)  under  the  condition  that  they  will 

be    finite   at   u   =    0    are: 

a  *°  N 


H    =    21       ^S^  T^  (o^u)   ®M   /  C22) 


rv\-=o 

and 


1  tvx-O 


A 


2-  <v^0  d^ 


H     "   2      B3ml«i<l*|u.)®«    7  (23) 


25 


Since  J  Cix)  =  i   I  Cx)  only  Eq.  C22)  can  be  used  without 
m  m 

losing  generality. 

Boundary  conditions  for  a  clamped  plate  are 

ur(\e)=  o 

Combining  EQ.  C2Q)  and  (22)  yields 


(24) 


A 


DO 


fc-^^+M^  $Aum)  u  w 


urM)=2  J^^U^) B^  ©„  (8-) 


C25) 


n\=.o 


where 


Applying  the  boundary  conditions,  Eqs .  C24)  to  (25)  and  the  re- 
quirement that  solutions  will  be  non-trivial  the  characteristic 
equations  are  obtained 


/ 


y& 


=  0 


(26) 


IL*1 


where  m  =  0,1,2,  ...  and  derivatives  are  w.r.t.  u.   The  dimension- 
less  frequencies-^-.,  m  =  0,1,2,  ...;  j  =  1,  2,  ...,  for  which 
Eq.  (26)  is  satisfied  are  the  eigenf requencies  to  which  there 


are  corresponding  values 


lmj  '   2mj '    3mj '   lm j 
corresponding  eigenfunctions  are 


and  *J  ~  .  .   The 
2m  j 
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A 

UJ". 
mJ 


Wmj  (0  ©M  (&) 


(t  r*  ■  >  jm(°W^t  (*  .  a)}©  [&) 


(a) 


*J 


m      —    0 /    1/    2,     ...;    j    —    1 /    2,     ...     . 


HU^,^)®^ 


-)<v0 


SJmk*ja>^«j£*(Su 


(V  -Q  Ja    '"^ 


x 


Ml    („/     .  ul)-^-    .   X,,,  (tf    .a) 


(27) 


(b) 


r  ,)i^^(r  _l)|^w\  ^«0 


"1  I^C^wj) -  ^  I^l^j) 


ill  U   /         -j-  /         ■>—   /         •••         /  I      ~™*       -L   /         — *    f  •••  * 

In  the  case  of  axisymmetric  motion,  m  =  0,  Eq.  (26)  takes 
the  form  (eliminating  the  derivatives) 


cLlCc*) 


c^t-o  *<  I  ct)  i;  &o  -(<^~!>  j.oo 


=    0  (28) 


and   the   eigenfunctions 
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A 


.IK-)-^}IoK^ 


v^"*  '"j^jr0^ 


IrM^.M^Mjt^'V^ 


\ 


=   0  (29) 

/ 
because   ^2/^9/=  0.   Thus,  there  is  no  coupling  between  U/~ 

or  QJ   and  IL)    in  ax i symmetric  motion. 

There  can  be  rotations  MJ  ^  0  that  are  axi symmetric,  but 

A  ^ 

they  will  be  independent  of  UT  and  11)  . 

The  characteristic  equation  for  independent  axi symmetric 

A 

(i)/v  motion  is  obtained  from  Eq.  (28)  by  setting  the  factor 
<*3  J,  C*$)  -  0.  (30) 

When  the  roots  of  Eq .  (30)  are  inserted  into  ^  g   the  _TL 

A 

values  are  obtained  for  axisymmetric  U/  motion.   The  ap- 
lied  force  is  normal  to  the  plate  here  and  these  modes  are  of 
no  interest.   This  can  be  seen  easily  after  calculating  a  value 

For  a  thick  steel  plate,  for  example,  —  =  0.2, XL  /   37  and 

a  r 

for  a  thin  plate  il  would  be  even  greater. 

Finally,  Eq.  (28)  for  this  problem  takes  the  form 

C  ^7-  1)*,  J0(  <XZ  )  Jx(^3  -  (  <£-  1^l  A,  (*.  )Ji(°(1)  -  0  (31) 

By  using  Eq.  (21a)  the  roots  -fl-  of  this  equation  can  be  ob- 

h  h 

tained.   Symmetric  eig  enf  requencies  for  —  =  0.04  and  —  =  0.2 

a  a 

are  given  in  Table  I  of  Ref.  16. 

ocVo  -n-V  '^riv/ 

N  when  \  (c°/c    )  C32) 
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Then   calling 

ft      =  —  o(      when  ^2    \^'    t^ie   characteristic    E<3 
ax i symmetric  motion   becomes 

and  corresponding  eigenfunctions  are 


(31)    for 


=    0 


(33) 


A 

or. 


J 


=   Jlfcxf 


c  «,;«■)  - 


0  v      I 


'J 


1.0* 


Ja 


A 


-c<r--i  jog  J,  (^aj-crj-uJe^'j)^  1,(6-11.)  (34) 

J  Xo^r) 


D.   COUPLING  OF  ACOUSTIC  FLUID  AND  PLATE 

The  following  developments  are  also  taken  from  Alper  and 
Magrab  [6].   On  the  5=0  surface,  which  represents  the  plate, 
boundary  conditions  are:   The  normal  fluid  velocity  and  the 
acoustic  pressure  on  the  fluid  are  equal  to  the  transverse 
plate  velocity  and  the  pressure  on  the  boundary  of  the  plate, 
respectively. 

Using  Eq.  (30)  , 


(35) 


1"° 


(36) 


The  external  force  (applied  force  per  unit  area  plus 

acoustic  pressure  at  the  plate  boundary  of  fluid)  is  given  as 
A 


(37) 
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Hence,  if  Eq.  (37)  is  substituted  into  Eq.  (19c) ,  then,  use  of 
Eqs .  (19)  and  (3  6)  couples  the  motion  of  the  plate  and  the 
acoustic  fluid,  so  that  solutions  for  yJ    and  VXT  can  be  found 
under  the  specified  force  distribution   1   .   When  Q  =  0  and 
SL   =  IL  .  ,  the  eigenfunctions   ( \jU".  ,  uj~   .  ,  llL    )  satisfying 
combined  Eq .  (19)  and  Eq .  (37)  could  be  found. 

Let  the  solution  for  the  plate  motion  be  expressed  as  expan- 
sions in  its  eigenfunctions : 
A     °° ,  0°  A 

ur=  zZ  Z.  D_uj;t 

A      **      g  ,  A 

<v\-0   J-\ 

Putting  Eq.  (3  3)  into  the  linear  and  homogenous  combination 
of  Eq.  (19)  and  (37)  there  is  obtained: 


The  orthogonality  condition  for  a  clamped  circular  Mindlin- 
Timoshenko  plate  is  defined  as 

=  ^mn  «"kj  0  "  ^o  J  n  ^mj  (40) 
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where   £)j_  is  the  Kronecker  delta  and 


\ 


N 


IT13 


•1  w^U^]-* 


(41) 


Using  Eq.  (27a)  and  performing  appropriate  integrations 


"M      *"J 


Z 


/ f  w  C  X2- 


^  ^  rJ  1 


-^ 


tf^ 


Nmr^G^, 


A 

OfAT 


VJ  =  ~*  V"  WJ '  V^,    -     .,>•     Hv,  fo*T  ''O  +  ^ *  G.-  ^ 


1(V\T  +"  I  *i*v>j 


4$ 


A  ".I 


6*v  k^-L^^sO-fL  ML^o 


M  >2MJ     '^   '  J/        Smj 


rv\ 


4- 


-V 


z  2.      "^+i^^j'%j; 


^i<v\j    +■  l^wnj 


-for  Si}  ( 


T 


Co  ^ 
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in  which 


tx=uh£ 


Z\JBL 


«v\j         " v) 


2J*J 


A 


*J 


<*  Tm  C*^)  -  *3nj  JM-H  (?<&mj) 


SMW-  im1c-)-jm.,(-)Jm+i^) 


J  i 


G(.Cx)^[iM1(x)-^u)lM.iCixo]i: 

—  - 


A 


H  „,(.■*-.,  *z.")  = 


*,T«CMJMtl(*0^lTMooy»*d) 


z 


For  the  axisymmetric  case,  m  =  0 
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A 

Multiplying    the    first   of    Eq.     (39)    by    X"      Tuivif11^®'    ' 

the   second   by   ^    ^     .  IX  dud©-  t    and   the   third   by 

A  '    J 

\xT-    U.  cLw.    a  9    t    adding   the   three   equations   and    integrating   over 

0  ^  Q  ^  ?_jt      f  04  W-vM      and  using   the   orthogonality   condition 

Eq.     (40)     leads    to 


Iff  J 


J 


\T 


1         A 

quT.  U-duiQ    (4  2 

1      iv>T 


0    o 


m  =    0,1,2,...;    j    =    1,2 

By  expressing  f(u,^)  in  a  Fourier  series  f(u,T7)  = 


(v\=Q 


Fm(u)  = 

m 


in 


0+O* 


^9-j  <3)M  O)  ^ 


(43 


and  using  Eq.  (12)  to  express  \U   there  results 


Combining  Eq.  (4  2) ,  (43)  and  (44)  and 
with  respect  to  w      yields 

VC   /        <*-i     wit  ^  «t/  — 


(44) 
performing    integration 


J  ^  «*,rtt 

m  =    0,1,2,.,.;    t   =    1,2,     ...     . 


(45) 
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»h*~       Vj  =/  \X/MJU)  ^  (_iil^)  a  du 


and 


"wj  -j    ^m  U)  \X4f  CO  *  .  All 


Expressions  for  g     are 


d<*,nj 


V 


9      *N   Q  2 


where 


3 


60         Co 


a. 


<H> 


1 


(46) 


(47) 


V  2- 


$Z  =   1,2  and 


Pc 


Pd 


|  N*1    MO  ^ 


l«i*l 


J    £o   ^  ^|!(-W|)!(k+*i41+i)\   l 


J 


2a+M 


For  the  axisymmetric  case  m  =  0. 

Combining  Eq.  (12)  and  (3  6)  and  using  the  first  of  Eq.  (38) 
for  w  yields 


The  functions   S^  '  (  -  iJL,V  )  ,     (n-m)  even,  are  orthogonal 
on  the  interval  0  4^-^1  *   The  orthogonality  condition  is 
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I-0 

(48) 


(49) 


(n-m)  even;  m  =  0,1,2,  ...  . 


where 


(50) 
(n-m)  even. 

Multiplying  Eq.  (48)  by  S^~'  ,  (  -  ill,  01  )  ^a1"^,  integrating  over 
the  interval  0  4^i  I  '  solving  for  A  and  using  the  definition 
Eq.  (46)  yields: 


fc> 


^Snn  " 


t^  ^'c-i^o)^,  «j<"v\r 


7  i  -^t  9 


(51) 


J=l 


m   =  0,1,2,  ...  ;  n=m,  m+2,  m  +  4. 
Putting  Eq.  (51)  into  Eq .  (45) 


52) 


^  r?-\L 


This  represents  an  infinite  number  of  sets  of  infinitely 
long  linear  non-homogenous  equations  in  the  coefficients  D  • 
For  each  value  of  m,  there  is  a  set  of  equations,  t  =  1,2, 
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The  infinite  set  of  equations,  yields  solutions  for  D  .;  i  =  1, 

2,  ...  which,  when  put  Eq.  (38)  gives  the  displacements  of  the 

plate  and,  when  put  into  Eq.  (51),  gives  the  coefficients  A 

The  coefficients  A   are  then  used  to  determine  the  m   harmonic 

mn 

of  the  potential  field  using  Eq.  (12)  and  thereby  the  m   har- 
monic of  the  pressure  using  Eq.  (5) . 

Instead  of  using  an  infinite  set  of  equations  there  is  a 
necessity  to  use  truncation  with  a  resulting  limitation  on 
accuracy  (N  =  number  of  terms  in  the  truncated  solutions) .  Alper 
[16]  found  that  for  the  far  field  pressure  using  N  =  11  for  thin 
C|-  0.04)  and  N=  16  for  thick  (|  .  Q_2)  plates  gives  solutions 

which  are  accurate  to  within  0.2  percent. 

Alper  calculated  patterns  for  the  first  three  resonant  fre- 

quencies  for  thin  (—  =  0.04)  and  thick  (—  =  0.2)  steel  plates. 

a  a 

These  are  presented  in  Figures  3  and  4 .   Calculated  pressures 
in  the  far  field  on  the  plate  axis  as  a  function  of  dimension- 
less  frequency  for  thin  and  thick  plates  are  presented  in  Figures 
5  and  6.   Physical  properties  for  steel  plates  with  water  as  the 

fluid  medium   are:        k    =    0.86,    V  =    0.3,      —=    0.42,    ijn   =    0.13  2 

T  £-5  J 

and  P0  =  1 . 

As  a  result  it  was  decided  that  the  best  choice  for  building 
a  transducer  to  meet  the  desired  specifications  would  be  a  thick 
plate.   A  sketch  illustrating  the  form  of  plate  displacement  for 
the  first  (Sl=   1.9)  and  second  (-0-=  6.9)  symmetric  modes  is 
presented  in  Figure  7 . 

The  computer  programs  used  for  the  theoretical  calculations, 
which  were  designed  for  a  CDC  64  00  computer,  were  obtained  from 
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Dr.  Alper.   These  have  been  modified  for  the  CDC  6  50  0  and  are 
listed  in  Appendix  A.   A  shortage  of  time  prohibited  actual  use 
of  these  modified  programs  for  design  work. 
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III.   EXPERIMENTAL  PROCEDURES 

An  experimental  test  of  the  calculations  of  Alper  and 
Magrab  was  conducted.   A  small,  clamped  edge  plate  transducer 
was  constructed  with  dimensions  which  corresponded  to  the  pub- 
lished calculated  results  for  the  thick  plate,  as  shown  in  Fig- 
ure 8.   The  material  used  is  steel.   The  ratio  of  plate  thick- 
ness to  radius,  h/a,  was  chosen  to  be  0.2.   The  thickness  was 
chosen  so  that  the  second,  in-fluid  resonance  mode  would  be  about 
75kHz.   This  corresponds  to  a  value  of  the  parameter  -H.  =  6.9. 
Exterior  dimensions  were  determined  by  the  positioning  and 
mounting  hole  on  the  cylindrically-shaped  torpedo  extension  sec- 
tion, which  served  as  the  acoustic  baffle  for  the  transducer. 
The  height  of  the  outer  cylindrical  section  at  the  rear  was 
chosen  to  be  about  3/4  wave  length  in  order  to  create  a  high 
acoustic  impedance  for  the  plate  edge  in  the  neighborhood  of 
75kHz. 

The  driver  for  this  transducer  consisted  of  a  piezoelectric 
ceramic  longitudinal  vibrator,  similar  to  the  sketch  shown  in 
Figure  9.   The  ceramic,  a  lead-zirconate  titanate  type,  con- 
sisted of  two  axially  polarized  circular  discs  with  a  hole  in  the 
center.   Masses  on  the  ends  and  a  clamping  screw  were  constructed 
and  their  dimensions  were  selected  to  make  the  half -wavelength 
(free-free)  longitudinal  resonance  to  occur  near  7  5kHz.   A 
silicone  grease  was  used  between  joints  to  ensure  good  acoustic 
coupling. 
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In  order  to  provide  a  more  nearly  point  source  of  excitation, 
the  area  of  contact  of  the  driver  with  the  plate  was  reduced 
by  the  use  of  a  small  nut  on  the  clamping  screw  between  the 
driver  and  the  plate.   Later,  measurements  were  made  with  the 
full  face  of  the  vibrator  clamped  to  the  plate. 

Resonance  frequencies  of  the  completed  transducer  assembly 
were  determined  using  measurements  of  the  electrical  admittance 
of  the  piezoelectric  elements.   A  Dranetz  Impedance-Admittance 
meter  was  used  for  this  purpose.   Typical  graphical  results  of 
conductance  versus  frequency  are  shown  in  Figures  12  through  14. 

The  construction  of  the  longitudinal  vibrator  used  in  the 
acoustic  measurements  prohibited  the  obtaining  of  a  satisfactory 
isolation  of  the  electrodes  from  electrical  ground  of  the  system, 
needed  for  making  electrical  admittance  measurements  with  the 
unit  in  water.   Therefore,  resonance  frequencies  were  not  mea- 
sured in  water.   However,  measurements  were  conducted  of  reson- 
ances of  the  driver  with  the  transducer  in  air  and  in  water 
using  a  different,  isolated  driver.   This  driver  was  a  radially 
polarized  cylinder  of  barium  titanate,  2.54  cm  long  by  1.27  cm 
outside  diameter.   It  was  clamped  to  the  plate  using  a  screw 
and  an  end  mass.   Although  the  free-free  half  wave  length 
resonance  frequency  for  this  element  differed  by  several  kilo- 
hertz  from  that  of  the  other  driver,  the  observed  resonances 
of  the  combined  systems  differed  only  slightly.   The  change  in 
resonance  frequency  caused  by  the  immersion  in  water  was 
within  the  uncertainty  of  frequency  setting.   The  effects  of 
radiation  damping  on  the  sharpness  resonance  was  noted. 
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The  acoustical  measurements  consisted  of  measuring  the  far 
field  acoustic  pressure  as  a  function  of  angle.   The  arrange- 
ment of  apparatus  is  shown  in  schematic  form  in  Figures  10-11. 
A  pulse  method  was  used.   Although  the  walls  and  floor  of  the 
water  tanks  are  essentially  anechoic,  the  reflection  from  the 
water-air  surface  can  be  significant,  necessitating  the  use  of 
a  box-car  integrator.   The  gate  of  the  box-car  integrator  was 
adjusted  so  that  only  the  sound  wave  traveling  by  a  direct 
path  from  the  source  to  the  hydrophone  was  measured. 

The  cylindrical  baffle  consists  of  an  aluminum  alloy  cast- 
ing about  2  6.25  cm  long  and  3  2.4  cm  outside  diameter.   Plastic 
end  caps  were  used.   The  wall  thickness  is  not  uniform  but  is 
about  2.2  cm.   A  sketch  is  shown  in  Figure  12. 

The  transducer  is  bolted  into  a  hole  which  is  located  as 
shown.   Lead  weights  were  added  inside  the  cylinder  to  make 
the  assembly  negatively  buoyant. 

The  rotation  system  could  be  arranged  to  measure  the  radia- 
tion pattern  about  either  the  longitudinal  axis  of  the  cylindej 
or  the  perpendicular  axis.   Measurements  about  both  axes  were 
made.   The  mechanical  drive  used  to  rotate  the  transducer  in- 
cludes a  potentiometer  on  the  gear  train  which  produces  a  D.C. 
voltage  proportional  to  angle  of  rotation.   This  voltage  is 
applied  to  the  x  axis  of  x-y  plotter.   The  calibrated  LC-10 
hydrophone  used  as  a  receiver  and  the  transducer  were 
separated  by  a  distance  of  1  m  with  the  hydrophone  and  the 
center  of  the  transducer -at  a  depth  of  1.92  .   The 
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output  of  the  hydrophone  was  amplified,  rectified  by  the  envelope 
detector  and  fed  to  the  box-car  integrator.   In  order  to  elimi- 
nate the  effects  of  surface  reflections,  and  minimize  the  ef- 
fects of  standing  waves  in  the  tank  a  pulsed  signal  was  used.   The 
box-car  integrator,  gated  to  the  direct  path  pulse  produces  a 
D.C.  voltage  that  is  proportional  to  the  average  amplitude  of 
the  received  pulse  envelope.   This  D.C.  voltage  is  fed  to  the 
logarithmic  converter  which  generates  the  ordinate  of  the  beam 
pattern.   The  attenuator  is  used  to  calibrate  the  y  axis  to 
provide  a  decibel  scale.   It  is  needed  because  of  the  non- 
linearity  of  the  envelope  detector. 

Each  radiation  pattern  was  calibrated  by  adjusting  the  x-y 
recorder  and  driving  voltage  to  obtain  at  least  half  scale 
reading  on  the  integrator,  and  then  by  the  means  of  attenuator, 
decibel  increments  on  the  y  scale  were  marked. 

Results  of  radiation  patterns  about  the  perpendicular  axis 
and  the  cylindrical  axis  are  shown  in  Figures  20  and  21  and 
Figures  22  and  23,  respectively.   These  were  made  with  a  small 
area  of  coupling  between  the  driver  and  the  plate.  Then,  with 
the  small  brass  nut  removed  from  the  coupling  between  driver 
and  plate  the  area  of  contact  was  increased  to  that  of  the 
disk  at  the  end  of  the  driver.   Rotation  was  about  the  vertical 
axis-     Results  are  shown  in  Figures  24  through  26. 

The  box-car  integrator  gate  width  was  suitably  maintained 
to  ensure  that  the  integration  would  be  done  at  steady  state 
region  of  received  signal. 

After  some  experimental  work  was  done  it  was  noticed  that 
bubbles  produced  in  the  water  greatly  affect  the  results  and 
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reduce  loading  effects  of  the  fluid.   Hence  to  prevent  this 
effect  a  long  brush  was  used  for  both  hydrophone  and  transducer. 
Before  each  pattern  was  recorded,  the  source  and  receiver  were 
brushed  off  in  order  that  no  effects  due  to  bubbles  would  occur. 

A  second  series  of  measurements  using  the  composite  driver 
with  small  coupling  area  was  conducted.   The  objective  was  two- 
fold.  The  first  was  to  identify  the  modes  of  vibration  associ- 
ated with  the  several  resonance  frequencies  noted  previously 
and,  particularly,  to  determine  which  were  the  first  and  second 
modes.   The  second  objective  was  to  recheck  the  results  of  earli- 
er measurements. 

For  the  driver  alone  (free-free)  conductance  plotted  as  a 
function  of  frequency  is  shown  in  Fig.  13.   Two  resonances  are 
noted.   It  seemed  clear  from  observing  the  sensitivity  to  clamping 
provided  by  grasping  it  with  one's  fingers  that  the  uppermost 
resonance  involved  flexural  plate  resonance  of  the  lower  and 
mass.   The  lower  resonance  appeared  to  be  the  usual  free-free 
mode  of  vibration. 

When  this  driver  is  mounted  on  the  transducer,  a  number 
of  resonances  can  be  observed.   A  plot  of  driver  conductance 
as  a  function  of  frequency  is  presented  in  Figures  14  and  15, 
for  the  transducer  alone  and  for  the  case  where  the  transducer 
is  mounted  in  the  cylindrical  baffle. 

A  small  capacitive  displacement  probe  was  used  to  measure 
the  relative  displacement  amplitude  at  various  positions  on 
the  outer  face  of  the  transducer.   The  probe  consisted  essen- 
tially of  an  8  mm.  diameter  flat  electrode  which  could  be 
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positioned  very  close  to  the  plate  (about  0.1  mm) .   This  elec- 
trode and  the  plate  itself  formed  an  electrical  capacitor.   The 
electrical  circuit  provided  polarization  potential  for  this 
capacitor  and  measured  the  amplitude  of  the  alternating  poten- 
tial caused  by  changing  capacitance  when  the  plate  vibrated. 
Although  the  probe  measures  the  average  effect  of  change  in 
capacitance  over  a  distance  of  6  or  8  mm,  it  was  sufficiently 
sensitive  to  get  a  reasonable  estimate  of  the  relative  amplitude 
of  plate  motion  as  a  function  of  distance  from  the  center  of  the 
disk.   A  schematic  representation  of  the  displacement  measure- 
ment apparatus  is  shown  in  Figure  16. 

A  plot  of  the  relative  amplitudes  is  presented  in  Figures 
17  and  18  for  each  of  the  transducer  resonances  noted  previously 
The  motion  for  the  01  and  02  modes  was  re-confirmed  after  the 
transducer  was  clamped  in  the  cylindrical  baffle. 
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IV.   DISCUSSION  OF  RESULTS 

A.  EFFECTS  OF  AXIS  OF  ROTATION  ON  PATTERN 

Radiation  patterns  were  obtained  for  the  same  transducer 
but  with  different  axis  of  rotation,  one  parallel  to  the  cylin- 
der axis  and  the  other  perpendicular  to  the  cylinder  axis.   The 
results  for  two  such  measurements  are  presented  in  Figures  20- 
21  and  22-23. 

It  is  concluded  that  although  some  differences  between 
experimental  patterns  exist,  these  do  not  obscure  the  basic 
features  of  the  radiation  pattern  to  be  used  for  comparison 
between  theory  and  experiment.   For  purposes  of  convenience, 
all  remaining  patterns  were  made  by  rotating  the  baffle  about 
its  perpendicular  axis. 

B.  EFFECTS  OF  VARYING  AREA  OF  CONTACT  BETWEEN  DRIVER  AND  PLATE 
Some  radiation  pattern  measurements  were  made  in  which  the 

size  of  the  contact  area  between  the  longitudinal  vibrates  was 
increased  from  about  5  mm  diameter  to  the  diameter  of  the  driver 
end  mass,  17  mm.   Since  the  diameter  of  the  driver  end  mass  is 
comparable  to  the  dimension  of  the  first  nodal  circle  for  the  02 
mode  of  the  plate,  it  should  be  expected  that  the  details  of 
the  displacement  of  the  plate  at  the  higher  frequencies  would 
not  be  a  good  approximation  to  that  assumed  for  the  point  drive 
case.   The  radiation  patterns  made  at  high  frequencies  using 
this  larger  area  of  coupling  were  very  complicated  and  difficult 
to  interpret  and  they  are  not  included  here.   The  pattern  at 
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low  frequencies  made  with  the  larger  area  of  coupling  were  not 
significantly  different  from  those  using  the  small  coupling 
area.   These  are  presented  in  Figures  24  through  26. 

C.  IDENTIFICATION  OF  MODES  OF  VIBRATION  OF  THE  PLATE 

The  existence  of  a  number  of  resonances  in  the  neighbor- 
hood of  frequencies  calculated  for  the  01  and  02  modes  of  the 
transducer  plate  appear  to  be  due  to  the  frequency  dependent 
impedance  of  the  various  parts  attached  to  the  plate.   The 
driver  element  is  one  such  attachment.   However,  the  boundary 
conditions  at  the  plate  edge  are  probably  more  important.   The 
massive  cylinder  at  the  rear  of  the  plate  was  intended  to 
provide  a  high  mechanical  impedance  at  the  plate  edge  to  simu- 
late the  clamped  edge  condition.   An  outer  flange  of  thickness 
comparable  to  (but  larger  than)  that  of  the  plate  itself  was 
needed  to  secure  the  transducer  in  the  baffle.   There  were 
several  modes  of  vibration  which  involved  significant  vibration 
of  the  outer  flange.   Some  of  these  modes  were  damped  where 
the  transducer  was  clamped  into  the  baffle.   It  is  concluded 
that  the  resonances  at  31  and  63.9  kHz  were  associated  with 
the  clamped  edge  plate  modes  having  mode  number  01  and  02, 
respectively. 

D.  RADIATION  PATTERNS 

The  comparison  between  the  radiation  patterns  calculated 
by  Alper  [16]  and  those  measured  in  these  experiments  must  take 
into  account  several  factors  which  can  give  rise  to  differences 

First,  the  assumed  clamped  edge  boundary  condition  is  only 
approximated  in  the  tests. 
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Second,  the  baffle  used  is  cylindrical  in  shape  and  is 
therefore  neither  infinite  nor  plane,  so  the  assumption  of  the 
theory  of  an  infinite  plane  baffle  is  not  well  approximated 
even  though  the  diameter  of  the  cylinder  is  about  7.5  times 
the  diameter  of  the  plate. 

Third,  the  fact  that  the  baffle  is  not  perfectly  rigid  per- 
mits some  coupling  of  flexural  wave  energy  from  the  plate  into 
the  baffle  and  the  possibility  of  additional  radiation  of  sound 
into  the  water  from  the  baffle  itself. 

Fourth,  the  area  of  contact  between  the  driver  and  the  center 
of  the  plate  was  made  as  small  as  possible  but  it  was  not  a 
point  driving  force,  as  assumed  in  the  theory. 

The  comparison  between  theoretical  and  experimental  radia- 
tion patterns  should  best  be  made  by  inspecting  Figures  3  and  4 
and  27  and  29.   The  theoretical  patterns  due  to  Alper  (shown  in 
Figure  19)  have  been  redrawn  in  the  rectangular  log-linear  plot 
format  used  for  the  experimental  data- 

The  experimental  pattern  which  should  correspond  most  closely 
to  the  01  mode  is  Figure  27  at  a  frequency  of  about  31  kHz. 

Those  which  should  correspond  most  closely  to  the  02  mode 
are  shown  in  Figures  20-22-29  at  a  frequency  of  63  kHz. 

For  the  31  kHz  pattern  it  is  seen  that  there  are  some  small 
oscillations  of  about  3  dB    near  the  center  of  a  single  major 
lobe.   The  level  of  this  lobe  at  the  90°  angle  is  down  about 
18  dB  below  the  axial  value.   This  contrasts  with  the  theoretical 
90°  value  of  about  -2  dB.   The  reasons  for  this  disagreement 
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are  not  known,  but  it  is  believed  that  an  important  contribu- 
tor is  the  finite  size  and  curvature  of  the  baffle. 

For  the  63  kHz  patterns,  there  are  some  features  which 
correspond  well  to  the  theoretical  pattern.   The  angular 
position  of  the  first  minimum  is  at  about  23  degrees,  very 
close  to  the  theoretical  value   of  22.5  degrees.   The  level  of 
the  measured  maximum  of  the  size  lobe  is  about  the  same  as  that 
of  the  major  lobe.   The  theoretical  curves  indicate  a  side 
lobe  maximum  level  of  -2.5  dB  located  at  48°  and  a  level  at 
90°  of  -5.6  dB.   The  experimentally  observed  level  at  90°  is 
about  -20  dB.   Again,  this  lower  level  near  90°  is  believed  to 
be  due  largely  to  the  finite  baffle  dimension.   The  additional 
structure  on  the  measured  patterns  is  believed  to  be  due  to 
radiation  of  sound  from  the  baffle. 

Illustrations  of  the  complexity  of  radiation  patterns  at 
other  higher  frequencies  are  shown  in  Figures  21-23-3  0.   The 
frequencies  here  correspond  to  resonances  in  which  there  are 
additional  displacements  of  the  transducer  parts  particularly 
the  outer  flange  (see  Figures  16-17) .   Very  likely  the  baffle 
itself  was  also  excited  into  some  flexural  modes  of  vibration. 

It  is  concluded  that  the  major  features  of  the  radiation 
patterns  measured  agree  with  those  of  the  Alper  theory  and 
that  the  differences  are  probably  due  to  boundary  conditions  of 
the  plate  which  did  not  meet  those  assumed  in  the  theory. 

E.   RECOMMENDATIONS 

Since  the  goal  for  this  effort  was  to  develop  the  Alper- 
Magrab  theory  as  a  design  tool  for  the  development  of  a  flexural 
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disk  transducer  that  might  be  suitable  for  use  as  a  sound 
source  on  underwater  test  range  tracking,  the  following  comments 
and  recommendations  seem  to  be  appropriate. 

The  adaptation  of  Alper's  program  to  fit  the  CDC  6500  Com- 
puter, HAL,  is  virtually  complete.   Shortage  of  time  precluded 
completion  of  this  part  of  the  work.   It  is  recommended  that  the 
work  be  finished  and  that  it  be  used  as  intended. 

It  appears  that  some  sort  of  variation  on  the  radiation 
pattern  from  the  second  radially   symmetric  plate  mode  does  offer 
the  potential  for  controlling  the  source  level  near  the  axis 
and  yet  maintain  a  relatively  broad  beam  width  for  the  secondary 
lobes.   The  patterns  observed  here  show  a  deep  but  narrow  null 
at  about  22  degrees  off  the  central  axis.   It  may  be  possible 
that  this  null  could  reduced  by  proper  distribution  of  excitation 
force  on  the  plate.   It  also  is  possible  that  the  effect  of 
this  narrow  null  on  range  tracking  is  acceptable.   It  is  recom- 
mended that  the  experience  of  this  work  be  employed  to  con- 
struct and  test  on  the  tracking  range  a  flexural  disk  trans- 
ducer which  operates  in  this  mode. 

The  smoothest  radiation  pattern  is  given  by  the  plate 
operating  in  the  first  symmetric  mode.   The  transducer  used  in 
these  tests  was  not  optimized  for  this  mode.   It  is  likely 
that    a  transducer  designed  to  operate  in  this  mode  at  7  5  kHz 
would  have  such  a  small  plate  diameter  that  a  problem  would 
exist  in  fitting  an  adequate  sized  driver  into  the  available 
space.   If  the  program  development  tool  becomes  workable,  cal- 
culations should  be  made  to  test  the  design  at  7  5kHz  to  deter- 
mine feasibility. 
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A  much  more  likely  application  for  a  flexural  disk 
source  operating  in  the  first  mode  appears  to  be  at  frequencies 
significantly  lower  than  75  kHz,  such  as  are  used  on  other  test 
ranges. 


49 


Fluid    Side 


O    J    ^-Q 


Rigid 

BRFFLE 


Vrcuum    Side 


Figure  1.   Geometry  of  Problem 
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Figure  2.   Oblate  Spheroidal  Coordinate  System 
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Figure   17.      Relative  Plate   Displacement  of   Transducer 
(Unmounted,    in  Air] 


66 


RELATIVE  PLATE  DISPLACEMENT 

63.9  kH* 

OZ  MODE 


2> 


RELATIVE  PLATE  DISPLACEMENT 


31    kH* 

01  MODE 


r  (cm) 


Figure   18.      Relative   Plate   Displacement   for    01   and    02 
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Figure  20.   Radiation  Pattern  (small  area,  perpendicular 
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Figure  22,   Radiation  Pattern  (small  area,  perpendicular  axis) 
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Figure  25.   Radiation  Pattern  (large  area,  perpendicular  axis) 
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Figure  26.   Radiation  Pattern  (large  area,  perpendicular  axis) 
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Figure  28*   Radiation  Pattern  (small  area,  perpendicular  axis) 
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Figure  29.   Radiation  Pattern  (small  area,,  perpendicular  axis) 
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Figure  30.   Radiation  Pattern  (small  area,  perpendicular  axis) 
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APPENDIX  A 

COMPUTER  PROGRAMS  FOR  FAR  FIELD  PRESSURE 
AND  ASSOCIATED  PARAMETERS 

Alper ' s  original  programs  were  not  changed  very  much. 
For  the  main  program,  the  outputs  of  the  Program,  Roots,  and 
the  Program,  Caller,  (for  getting  sum')  are  needed.   After 
getting  the  output  of  Program,  Roots,  Program,  Caller  output 
is  obtained  by  using  these  results!   Finally,  the  main  pro- 
gram will  be  run  by  using  the  obtained  parameters  and  sums. 
For  this  reason  some  Formats  were  changed  to  fit  inputs  and 
outputs  each  other.   The  only  problem  which  still  was  not  solved 
is  that  of  Subroutine  Bessel. 
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